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Oncogenes such as K-ras mediate cellular and metabolic transformation during tumorigenesis. 
To analyze K-Ras-dependent metabolic alterations, we employed ^^C metabolic flux analysis (MFA), 
non-targeted tracer fate detection (NTFD) of ^^N-labeled glutamine, and transcriptomic profiling in 
mouse fibroblast and human carcinoma cell lines. Stable isotope-labeled glucose and glutamine 
tracers and computational determination of intracellular fluxes indicated that cells expressing 
oncogenic K-Ras exhibited enhanced glycolytic activity, decreased oxidative flux through the 
tricarboxylic acid (TCA) cycle, and increased utilization of glutamine for anabolic synthesis. 
Surprisingly, a non-canonical labeling of TCA cycle-associated metabolites was detected in both 
transformed cell lines. Transcriptional profiling detected elevated expression of several genes 
associated with glycolysis, glutamine metabolism, and nucleotide biosynthesis upon transforma- 
tion with oncogenic K-Ras. Chemical perturbation of enzymes along these pathways further 
supports the decoupling of glycolysis and TCA metabolism, with glutamine supplying increased 
carbon to drive the TCA cycle. These results provide evidence for a role of oncogenic K-Ras in the 
metabolic reprogramming of cancer cells. 
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Introduction 

Several decades ago Otto Warburg first described that tumors 
exhibit glycolytic metabolism with a reduced rate of oxidative 
phosphorylation, despite the availability of adequate oxygen. 
This phenomenon, known as the 'Warburg effect,' has been 
proposed to be a key driver of tumor progression (Warburg, 
1956). Since then, several researchers have observed altered 
glucose metabohsm associated with a high rate of lactate 
secretion in cancer cells and tumor tissues (Sauer et aU 1982; 
Steinberg et aU 1992). Furthermore, aerobic glycolysis 
represents a robust hallmark of cancer that is employed for 
tumor detection by positron emission tomography analysis 
with the 2-[18F]fluoro-2-deoxy-D-glucose tracer (Laing et aU 
2009). However, we and others have recently provided 
evidence on the important anabohc role of glutamine in tumor 
cell prohferation (DeBerardinis et aZ, 2007; Yuneva et aZ, 2007; 
Gagho et al, 2009). Indeed, by using different experimental 
approaches, it has been shown that glycolytic cancer cells 
consume more glutamine as compared with their normal 
counterparts to synthesize proteins, nucleotides and fatty 



acids and produce energy (Kovacevic and McGivan, 1983; 
Baggetto, 1992; DeBerardinis et aU 2007). 

Several reports suggest that both nutrient uptake changes and 
metabolic alterations are under direct control of ras or myc 
oncogenes (Chiaradonna et aU 2006b; Wise et aZ, 2008; Vander 
Heiden et aU 2009). In particular, oncogenic Ras proteins, 
identified in 25% of human cancers (Bos, 1989), correlate with 
metabolic alterations, including increased glucose and 
glutamine consumption, lactic acid accumulation, altered 
expression of mitochondrial genes, increased reactive oxygen 
species (ROS) production and reduced mitochondrial activity 
(Bos, 1989; Vizan etal, 2005; Chiaradonna etal, 2006a; Yun etal 
2009; Weinberg et al 2010). A result of this metabolic 
reprogramming is the dependence of K-Ras transformed cells 
on glucose and glutamine availability, since their withdrawal 
induces apoptosis and cell-cycle arrest, respectively (Rama- 
nathan et al 2005; Telang et al 2006; Yun et al 2009} . However, 
the precise metabolic effects downstream of oncogenic Ras 
signaling in cancer cells have not been completely elucidated. 

Due to the interconnected nature of metabolic pathways and 
the pleiotropic effects mediated by oncogenic signals, a 



© 201 1 EMBO and Macmillan Publishers Limited 



Molecular Systems Biology 201 1 1 



K-Ras decouples glucose and glutamine metabolism 
D Gaglio et al 



systems approach is required to elucidate the mechanisms of 
such transforming events. To date, most large-scale analyses of 
tumor cells have employed microarray technology that 
provides a robust means of analyzing transcriptional mod- 
ifications associated with physio/pathological changes (Ross 
et aU 2000; Scherf et aU 2000} . However, various levels of post- 
transcriptional controls may not be captured by these analyses 
(Mata et aU 2005; Metallo and Vander Heiden, 2010), and new 
approaches are now emerging to increase our knowledge 
about cancer cell physiology (Liotta and Petricoin, 2000; 
Laubenbacher et aU 2009; Kreeger and Lauffenburger, 2010). 

Metabolomic techniques offer a more direct means of 
studying metabolism at the systems level. Measuring metabo- 
lite concentrations, indeed, represents a sensitive approach that 
generates detailed 'snapshots' of biological processes (Hiller 
et aU 2009} . Nevertheless, metabolic routes are best character- 
ized by the measurement of fluxes, which describe the actual 
functionality of a given enzyme or pathway (Stephanopoulos 
and Vallino, 1991; Sauer, 2006). To this end, isotopic tracers and 
computational algorithms enable the quantitative estimation of 
intracellular fluxes and associated confidence intervals for a 
given system (Metallo et aZ, 2009; Hiller et aZ, 2010), and such 
methods are now effectively applied to mammalian cells (Vizan 
et al 2005; Munger et al 2008). 

Therefore, to better understand the regulation of cancer cell 
metabohsm and to identify key metabohc routes altered in 
K-Ras transformed cells, we have apphed a systems-level 
approach based on the combined apphcation of metabohc 
and transcriptional analyses. We have monitored the flux of 
^^C-labeled glucose and glutamine as well as [a^^N] glutamine 
into downstream metabolites in normal and transformed cells 
and performed a detailed comparison with the transcriptional 
profiles obtained from the same cell lines. Herein, we show 
that K-ms oncogene expression enhances glucose uptake but 
decreases its utilization in the tricarboxylic acid (TCA) cycle 
and associated anabohc pathways. Furthermore, we show that 
while K-Ras transformation decreases overall flux through the 
TCA cycle, it increases utilization of the carbon backbone and 
amino -nitrogen moieties of glutamine either through TCA 
cycle or transamination activities in order to sustain bio syn- 
thetic reactions, including amino-acid, nucleotide and glu- 
tathione synthesis. Finally, we present evidence describing the 
dependence of K-Ras on glutamine metabohsm, as inhibition 
of key enzymes along this pathway specifically compromises 
the prohferation of transformed cells. 

Results 

K-Ras transformation enhances glycolytic flux 
and decreases oxidative TCA cycle flux in 
mouse fibroblasts 

We have previously shown that transformed NIH3T3 mouse 
fibroblasts expressing an oncogenic K-Ras protein (G12V} 
(Shih etal 1981; Pulciani etal 1982; Bossu etal 2000} exhibit 
elevated sensitivity to glucose availability, reduced mitochon- 
drial function linked to a decrease of Complex I activity and 
protein level, and reduced prohferation in response to glucose, 
glutamine or galactose shortage as compared with immorta- 
lized NIH3T3 mouse fibroblasts (Chiaradonna et al 2005, 
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2006a; Gagho etal 2009; Baracca etal 2010). NIH3T3 cells are 
a genetically well-characterized immortalized cell line that has 
long been estabhshed as a model of 'normal' cells for the study 
of cell transformation, as these cells undergo contact inhibi- 
tion, exhibit no growth in soft agar, and do not form tumors in 
immunocompromised mice, in contrast to isogenic trans- 
formed lines (Bossu et al 2000). 

Therefore, we performed basic metabohc analyses to 
compare NIH3T3 cells (normal, N) and NIH3T3 K-Ras 
(transformed, T) 54 h (754) after plating. This time point 
was chosen since it corresponds to the time in which the two 
cell lines begin to show a different prohferative ability 
(Chiaradonna et al 2006b}. As shown in Figure lA, cell- 
specific glucose (Glc) and glutamine (Gin) consumption as 
well as lactate (Lac) and glutamate (Glu) secretion were 
measured. Significant increases in Glc uptake and Lac 
production in T cells were observed (Figure lA), consistent 
with other reports hnking ms mutations to aerobic glycolysis 
(Elstrom et al 2004; Vizan et al 2005). Although previous 
studies have shown an increase of Gin uptake in other kinds 
of cancer cell lines (Yuneva et al 2007; Wise et al 2008; 
Gao et al 2009), we did not detect any significant changes in 
Gin uptake and Glu secretion between the two cell lines. To 
further define N- and T-cell metabohsm, we have determined 
the relative abundance of various intracellular amino and 
organic acids by performing a targeted metabolomic analysis 
(Figure IB). Notably, metabolites associated with glycolysis 
and biosynthetic processes were present in larger quantities 
in T cehs, including pyruvate (Pyr), Lac, citrate (Cit), and 
aspartate (Asp) . On the other hand, TCA cycle intermediates 
such as a-ketoglutarate (AKG), succinate (Sue), and 
fumarate (Fum), which arise primarily from oxidative 
metabohsm, were less abundant in T cehs as compared 
with N cehs. 

Isotopic tracers provide an effective means of determining 
the relative contribution of a substrate to a particular 
metabolic pathway, and mass spectrometry enables quantifi- 
cation of the relative abundance of molecules with different 
numbers of isotopic label in each metabolite pool (i.e., MO, 
Ml, M2, ... where M2 abundance includes ah molecules with 
two ^^C-labeled atoms). For instance, uniformly labeled Glc 
([U-^^Cglglucose] is converted to M3 Pyr (with three ^^C 
atoms) and subsequently M2 Cit through the combined action 
of pyruvate dehydrogenase and citrate synthase. In the same 
way, M2 Cit, progressing along the TCA cycle, generates M2 
AKG and Glu by TCA cycle-related reactions (i.e., transamina- 
tion and reductive amination; Supplementary Figure SI). We 
cultured N and T cehs in the presence of ^^C-labeled glucose 
and quantified the mass isotopomer distribution (MID) 
of Cit and Glu via gas chromatography/mass spectrometry 
(GC/MS). These data indicated that T cehs oxidize less Glc 
carbon to Cit and Glu, presumably diverting more Pyr to Lac, 
as evidenced by the significant decreases in M2 labeling 
observed in these metabolites (Figure IC). 

Consistent with the observation of decreased Glc oxidation, 
parallel [U-^^Cs] glutamine labeling experiments indicated that 
Gin carbon utilization for metabohc intermediate synthesis 
(e.g. Asp) was increased in Tcehs. We specifically observed a 
relative increase in M3 and M4 Asp abundance (Figure ID), 
which are generated from M5 AKG downstream of reductive 
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Figure 1 K-Ras transformed fibroblasts decouple glycolysis and the TCA cycle. (A) Extracellular uptake and secretion of Glc, Lac, Gin, and Glu in N and T cells for 54 h 
of growth. (B) Relative metabolite abundances in N and T cells measured by GC/MS. (C) M2 labeling of Cit and Glu in N and T cells cultured with a 1:1 mixture of 
[1-^C]glucose and [U-''^C6]glucose. (D) Mass isotopomer distribution (MID) of Asp in cells cultured in the presence of [U-''^C5]glutamine. Error bars indicate s.e.m. 
(n=3). (E) Schematic representation of the metabolic routes involved in M3-M4 Asp labeling by using [U^^Cslglutamine tracer. In the scheme are represented the key 
cytoplasmic and mitochondrial enzymes involved in these metabolic routes. (F) Absolute expression values of some genes, involved in the pathways described in the 
scheme (E), in N and T cell lines. The two cell lines, grown in 25 mM Glc + 4 mM Gin (normal medium), were collected for transcriptional analysis at two different time 
points (48 and 72 h). The list of gene abbreviations is available in Supplementary information. 



carboxylation by isocitrate dehydrogenase and oxidation by 
oxoglutarate dehydrogenase, respectively (Supplementary 
Figure S2} . These results support the notion that Gin entering 
the TCA cycle of T cells may follow both the canonical forward 
oxidative reaction and a reductive carboxylation (Yoo et aU 
2008). Moreover, given the significant role of Asp in amino- 
acid and nucleotide synthesis, these results suggest that T cells 
are strongly dependent on Gin carbon for biosynthesis 
(Weinberg et aU 2010). 

Mass isotopomer data alone provides a somewhat limited 
view of cellular metabohsm. To better characterize N- and 
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T-cell metabohsm at the systems level, we integrated 
our extracellular uptake and secretion results with MID data 
obtained by culturing the cells with [U-^^Cs] glutamine. 
Using these data, we estimated intracellular fluxes in each 
cell type using an elementary metabolite unit (EMU) -based ^^C 
metabohc flux analysis (MFA) algorithm (Young et aU 2008). 
Upon obtaining acceptable fits, we performed a sensitivity 
analysis to generate 95% confidence intervals for each flux in 
the network (Antoniewicz et aU 2006; Antoniewicz et aZ, 
2007} . Flux results for N and T cells are Hsted in Table I, and the 
full dataset and model description are provided in Supple- 
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mentary information and Supplementary Figures S3 and S4. 
These data indicate that transformation induces a significant 
reorganization of central carbon metabohsm. T cells reduce 
flux through the PDH complex by ~50% (i.e., flux values of 
Pyr ^ AcCoAmit + CO2} . Furthermore, net fluxes throughout 
the TCA cycle were decreased by the same amount. Each of 
these changes was highly significant, as evidenced by the lack 
of overlap in confidence intervals between the two cell types. 
Although most exchange fluxes in reversible reactions of the 
TCA cycle could not be accurately determined, we observed no 
significant change in reductive flux through IDH or malic 
enzyme. Consistent changes in metabohc fluxes were 
observed when performing MFA using labeling data from cells 
cultured with ^^C-glucose (see Supplementary Information for 
detailed results), though the observability of some fluxes 



(e.g. IDH exchange flux) was significantly decreased when 
using these tracers. These findings provide evidence that 
transformation causes cells to decouple glycolysis and the TCA 
cycle, with an overall decrease in activity of oxidative reactions 
within the TCA cycle. 

Transcriptional profiling of N and T fibroblasts is 
largely coherent with metabolic flux analysis 

Various reports have indicated that oncogenes alter the 
expression of genes encoding metabohc enzymes (Mathupala 
et aU 1997 \ Zhong et al 1999). To identify oncogenic 
K-Ras-dependent transcriptional effects, we performed micro- 
array transcriptional profiling of N and T cells grown as above 
for 48 and 72 h. 



Table I Experimentally determined fluxes and 95% confidence intervals estimated from MFA model using [U-^^Cs] glutamine as a tracer 



Reaction 



Flux 



Normal 



Confidence interval 



Flux 



Transformed 



Confidence interval 



T/N 



Glycolysis 
GluCext^G6P 
G6P^F6P 
G6P^F6P 
F6P^DHAP + GAP 
DHAP^GAP 
DHAP^GAP 
GAP^3PG 
GAP^3PG 
3PG^Pyr 
Pyr -> Lac 
Pyr ^ Lac 
Lac — ^ LaCgxt 

Pentose phosphate pathway 
6PG->P5P + C02 
P5P + P5P^S7P + GAP 
S7P + GAP^F6P + E4P 
P5P + E4P^F6P + GAP 

TCA cycle 
Pyrmit + C02^0AA 
Mal^Pyrn,it + C02 
Mal^Pyrmit + C02 
Pyrmit^AcCoAn,it + C02 
AcCoAmit + OAA^Cit 
Cit^AKG + C02 
Cit^AKG + C02 
AKG^Suc + C02 
Sue Fum 
Sue Fum 
Fum Mai 
Fum <^ Mai 
Mal->OAA 
Mal^OAA 
Pyr->Pyrniit 
Pyr^Pyr^it 

Glutamine metabolism 
GlUgxt Gin 
Gln^Glu 
Glu GlUgxt 
Glu^AKG 
Glu ^ ARC 

Pyrmit + Glu ^ Ala + ARC 
OAA + Glu ^ Asp + ARC 



65.5 
62.5 
4.1 
64.2 
64.2 
36.7 
129.2 
0 

129.2 
114.9 
300.5 
114.9 



3.00 
0.83 
0.83 
0.83 



0 

1.8 
2.8 
16.1 
16.1 
10.1 
3.3 
12.1 
12.1 
37.2 
12.1 

0 
10.3 
36.2 
14 
74 



20 
20 
18 
2 

168 
0.15 
0.18 



61.2, 
58.2, 
0.0, 
59.9, 
59.9, 
0.0, 
120.6, 
0.0, 
120.6, 
105.9, 
0.0, 
105.9, 



69.9 

66.9 

Inf 

68.5 

68.5 

Inf 

137.9 

Inf 

137.9 

123.4 

Inf 

123.4 



2.71, 3.29 
0.73, 0.93 
0.73, 0.93 
0.73, 0.93 



0.0, 3.2 
0.9, 5.9 
0.0, 3.2 
13.5, 19.1 

13.5, 19.1 
9.1, 11.2 
2.7, 3.9 

10.6, 13.7 
10.6, 13.7 

0.0, Inf 
10.6, 13.7 

0.0, Inf 
6.7, 11.3 

0.0, Inf 
12.0, 17.0 
39.0, 133 



19.0, 22.0 
19.0, 21.0 
16.0, 19.0 
1.0, 3.0 
115, 341 
0.13, 0.16 
0.16, 0.20 



83.1 
82.1 
18.6 
82.4 
82.4 
13.1 
165 
0 
165 
158.1 
9.2 
158.1 



1.00 
0.17 
0.17 
0.17 



0 

1.8 
1.6 
8.6 
8.6 
5.4 
3.7 
7.4 
7.4 
0 

7.4 

42 250 

5.6 
57630 
7 

276 



16 
16 
14 
2 
202 
0.15 
0.18 



77.6, 
77.6, 
0.0, 
76.9, 
76.9, 
0.0, 
154.0, 
0.0, 
154.0, 
146.9, 
0.0, 
146.9, 



88.6 

87.6 

Inf 

88.0 

88.0 

Inf 

176.1 

Inf 

176.1 

169.2 

Inf 

169.2 



0.90, 1.10 
0.13, 0.20 
0.13, 0.20 
0.13, 0.20 



0.0, 1.9 
0.7, 4.7 
0.0, 1.9 
7.6, 9.6 
7.6, 9.6 
4.9, 6.0 
3.2,4.2 
6.1, 8.7 
6.1, 8.7 
0.0, 1.5 

6.1, 8.7 
1.5, Inf 
3.5, 6.1 

3.2, Inf 
6.0, 8.0 

87.0, 4566 



15.0, 17.0 
15.0, 17.0 
12.0, 15.0 
1.0, 3.0 
106, 1895 
0.14, 0.17 
0.16, 0.20 



U 

U 
ND 

U 

U 
ND 

U 

ND 

U 

U 
ND 

U 



D 
D 
D 
D 



NC 
NC 
NC 

D 

D 

D 

NC 
D 
D 

ND 
D 

ND 
D 

ND 
D 
U 



D 

D 

D 
NC 

U 
NC 
NC 



Detailed description, assumptions, complete results, and data are available in Supplementary information. Net flux: (uf-dr) (^) and Exchange flux: min (uf,i3r) (^). 
Some scaling and undetermined exchange fluxes (i.e., confidence intervals of infinity, Inf) are omitted for brevity. Comparison of fluxes (Transformed versus Normal) is 
described as U=up, D=down, NC=no change, and ND=no detectable. Fluxes are listed x 10^"* mol/(cell x h). 
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These data, represented by heat maps in Supplementary 
Figure S5, indicate that several genes of the investment phase 
of glycolysis {Hkl and Pfkp) as well as Pfkfb 2, 3 and 4, which 
through the regulation of fructose-2,6-biphosphate level 
allosterically control 6-phosphofructo-l -kinase, were consis- 
tently expressed at lower levels in T cells as compared with N 
cells (Supplementary Figure S5A} (green color). In addition, 
five genes encoding glycolytic enzymes [Gpil, PfkU Pfkm, 
AldoA and Tpil) were highly expressed, particularly at 48 h, in 
T cells as compared with N cells (red color) (Supplementary 
Figure S5A) . Moreover, a high expression for genes involved in 
the pay-off phase of glycolysis, where ATP is generated, was 
observed. Indeed, almost all the genes of this phase, at both 
time points 48 and 72 h, were expressed at a higher level in 
T cells as compared with N ones (Gapd/i, Pgkl^ Enol, Eno2, 
Pkm2, Ldha and Ldhb). 

Consistent with the reduction of TCA flux that we observed, 
three genes encoding kinase that inhibit the PDH enzyme 
complex activity [Pdkl, Pdk3, and Pdk4) were significantly 
expressed in T cells (Supplementary Figure S5B). In addition, 
expression of Ogd/i, which encodes one subunit of the 
oxoglutarate dehydrogenase complex, another key control 
point in the TCA cycle, was decreased in T cells (Supplemen- 
tary Figure S5B). On the contrary, few genes encoding for 
positive regulators of TCA cycle, as Aco2, SuclaZ and SucglZ, as 
well as Idh2, which is able to catalyze the reductive 
carboxylation of AKG, were expressed at higher levels in 
T cells throughout the time window of the analysis as 
compared with N ones (Supplementary Figure S5B). Taken 
together, these results indicate that numerous genes regulated 
by oncogenic K-Ras may drive the observed metabohc 
reprogramming, although we cannot exclude the possibility 
that other mechanisms (e.g., post-translational modifications, 
mass action, and allosteric enzyme regulation) also participate 
in the observed metabohc alterations of transformed cells. 



Transcriptional analysis of glutamine entry 
into TCA cycle 

Previous experiments (Figure 1) indicated that Gin contributes 
to precursors of cellular biomass (e.g., aspartate) to a greater 
extent in T cells as compared with N cells, despite the observed 
small decrease in uptake. These findings are consistent with 
other reports describing Myc-driven increases of Gin usage 
through glutaminolysis, mediated in part via expression of 
genes involved in Gin uptake and metabohsm (Wise et aZ, 
2008; Gao et aU 2009). To investigate this feature, we have 
analyzed the expression of genes encoding enzymes involved 
in Gin uptake and metabolism, with particular attention 
to those involved in metabohc pathways able to generate Asp. 
As shown in Figure IE, Asp can be produced from Gin both in 
cytoplasm and in mitochondria through AKG by oxidative 
decarboxylation and reductive carboxylation. Gin enters 
into the cells by high-affinity glutamine importer SlclaS and 
is metabolized to Glu by glutaminase [Gls) (Figure IE). 
Subsequently, Glu may be converted to AKG, either via 
transaminase activity of cytoplasmic and mitochondrial 
aspartate and alanine aminotransferases [Gotl and Got2 and 
Gptl and Gpt2, respectively) or via glutamate dehydrogenase- 
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dependent oxidative deamination [Gludl) (Figure IE). As 
shown in Figure IF, SlclaS, Glsl, and Gludl were consistently 
expressed at lower levels in T cells (about two-fold) relative to 
N cells. On the contrary, Gotl and 2, showed an opposite 
behavior, with both genes slightly upregulated in T cells 
compared with N cells. These transcriptional changes 
may partially justify the slight reduction in Gin uptake 
(e.g., downregulation of SlclaS in transformed cells) shown 
in Figure lA and suggest, by analysis of the metabohc routes 
depicted in Figure IE, a relevant role of AKG for the generation 
of labeled M4 Asp through oxidative metabohsm (see also 
transcriptional data in Figure IF and Supplementary Figure 
S5). In addition, in T cells, higher expression of Idh2 
(Supplementary Figure S5) and in particular of Ady, whose 
level was at least two-fold higher than N cells (Figure IF), is 
coherent with the increased M3 Asp labeling previously 
shown. However, while transcriptional data are for the most 
part consistent with the reduction of glutamine uptake and its 
use for biosynthesis, not all metabohc fluxes are expected to 
be only regulated at the transcriptional level. 



Transformed cells increase their use of glutamine 
nitrogen for biosynthetic processes 

To further investigate the role of Gin in cellular biosynthesis, 
we performed a non-targeted tracer fate detection (NTFD) 
analysis using amino-labeled glutamine ([a-^^N] glutamine) to 
identify and quantify metabolites containing glutamine- 
derived nitrogen. This method enables the detection and 
quantification of label in ah observable metabolites down- 
stream of the tracer substrate (Hiher et aU 2010) . Cehs require 
nitrogen atoms to build molecules such as nucleotides, amino 
acids, amino sugars with Gin and Glu representing the primary 
cellular reservoir of nitrogen. The amino nitrogen of Gin is 
retained on Glu and used as a nitrogen donor in transaminase 
reactions. As expected, the NTFD algorithm identified many 
amino acids and associated compounds that were highly 
labeled from [a-^^N] glutamine tracer, including Glu, Ala, Asp, 
and 5-oxoproline (5-oxo) (Figure 2A). With the exception of 
Asp, each of these metabolites was more significantly labeled 
in Tcehs as compared with N cehs. In addition, we observed 
that Ml and M2 labeling on the nucleobase adenine, generated 
from [a-^^N]glutamine-derived nitrogen atoms from Asp, was 
over 2-3 -fold higher in Tcehs compared with N control cehs 
(Figure 2B). 

Transcriptional data support results from NTFD analysis 
with [a-^^N] glutamine. Absolute gene expression values of N 
and T cehs and their ratio (T/N ratio) as weh, shown in 
Figure 2C, Supplementary Table SI, highlight the elevated 
expression of genes closely related to Gin utilization in Tcehs. 
In particular, the expression of genes involved in glutathione 
metabohsm (e.g., Gss, Gclm), hexosamine metabohsm 
(e.g., Gfptl), and nucleotide biosynthesis (e.g.. Cad, Adss, 
and Adsl) was increased in T cehs as compared with N cehs. 
Such transcriptional changes suggest increased activity of 
glutathione pathway and the enhancement of nucleotide 
synthesis, as indicated respectively by downstream metabolite 
5-0X0 labeling (Figure 2A) and Ml and M2 labeling on the 
nucleobase adenine (Figure 2B). In this regard, detailed 
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Figure 2 Labeled metabolites identified by NTFD analysis in cells cultured with [a-''^N]glutamine for 54 h. (A) M1 labeling on amino acids and derivative metabolites in 
N ( □ ) and T ( ■ ) cells. (B) M1 and M2 labeling observed in free adenine in N ( □ ) and T ( ■ ) cells. (C) Heat maps of genes involved in Gin metabolism of N (blue line) 
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analysis of genes specifically involved in purine metabolism 
(Figure 2D} have indicated a more sustained expression in T 
cells as compared with N cells. Altogether, these findings 
further suggest the main role of Gin in supporting anabohc 
processes of T cells as compared with N cells. 

K-Ras transformed human cells show metabolic 
alterations similar to those observed in mouse 
transformed cells 

To test whether human cancer cells exhibit similar metabolic 
alterations to the previously described mouse model of 
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transformation, we next sought to expand our metabolic and 
transcriptional analyses to the MDA-MB-231 human breast 
carcinoma cell line, which express an activated K-Ras protein 
(G13D}. To assess the glycolytic phenotype of MDA-MB-231 
cells, asynchronous cells grown in high and low Glc (25 and 
ImM, respectively), were followed in a time course of 96 h. 
As shown in Figure 3A, MDA-MB-231 cells grown in 25 mM Glc 
continued to prohferate and reached a much higher cell 
density compared with those in 1 mM Glc. In fact, in low Glc 
cells completely lost their prohferation ability and showed 
evident signs of cell death at 48 h, as revealed by prohferation 
curves (Figure 3A} and by morphological analysis (data not 
shown). In addition, MDA-MB-231 cells are dependent upon 
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Gin for proliferation; even low levels of Gin (0.5 mM) 
significantly limited growth compared with normal culture 
(Figure 3 A) . Since the effects of Glc and Gin shortage on cell 
viability and prohferation of MDA-MB-231 cells were very 
similar to those observed in the mouse model of K-Ras- 
dependent transformation, we analyzed the metabohc pheno- 
type of these human cancer cells in greater detail. Extracellular 
flux quantifications indicated that MDA-MB-231 cells convert 
Glc to Lac at approximately stoichiometric levels (i.e., 2:1; 
Figure 3B}. Indeed GC/MS analysis of cells cultured with a 1:1 
mixture of [l-^^C]glucose and [U-^^Celglucose have demon- 
strated that at least 90% of Lac was derived from Glc (i.e., M3 
Lac abundance was measured at 44 . 6 ± 0 . 1 % } . Compared with 
murine N and T cells, MDA-MB-231 cells consume much higher 
amounts of Glc relative to Gin (Figure 3C}. However, direct 
oxidative conversion of Gin still accounted for as much as 40 % 
of the Fum, Mai, and Asp pools in cells cultured with 
[U-^^Cs] glutamine (see M4 mass isotopomer abundances; 
Figure 3D}. Taken together, these data indicate that, much like 
transformed mouse cells, MDA-MB-231 cells divert most Glc to 
Lac and rely on Gin carbon for anaplerosis. 

Using extracellular flux data and MIDs obtained from MDA- 
MB-231 cells cultured for 54 h with [U-^^Cs] glutamine, we 
performed MFA to determine intracellular fluxes throughout 
central carbon metabohsm. The results of this analysis 
(Supplementary Figure S6; Supplementary Table S4} indicated 
that, in a similar way to mouse T cells, MDA-MB-231 cells 
exhibit high glycolytic flux and significantly lower TCA cycle 
activity. This comparison is best demonstrated by plotting the 
ratio of PDH flux to Glc consumption (Figure 3E}. While Pyr 
oxidation flux is ^25% of Glc uptake in control N cells, cells 
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with oncogenic K-Ras exhibit a PDH flux ^10% of Glc 
consumption. Therefore, oncogenic K-Ras consistently drives 
the decoupling of central carbon metabohsm, increasing the 
glycolytic flux to Lac and decreasing oxidative TCA flux in cells 
with different genetic backgrounds (Figures 1 and 3). 
Importantly, MDA-MB-231 cells carry a mutant p53 (Olivier 
et aU 2002), whereas K-Ras-NIH3T3 cells have wild-type p53 
(Wadhwa et aU 1999). Given the fact that recent literature 
results have imphcated the p53 tumor suppressor protein as a 
regulator of mitochondrial respiration and the Warburg effect 
(Matoba et aU 2006), we tested the effect of p53 status on 
these metabohc changes to ensure our results were specific to 
oncogenic K-Ras. Comparing p53^^^ wild-type (HCT116 
p53 + / + } or a p53"/~ KO (HCT116 p53~/"} colon carcinoma 
cells, which express oncogenic K-Ras (Bunz et aU 1998), we 
observed no significant differences in proliferation under 
nutrient deprivation (Supplementary Figure S7A and B) or 
intracellular metabolite abundances (Supplementary Figure 
S7C-F}. These data indicate that oncogenic K-Ras, at least in 
these cellular models, contributes to cancer-specific metabohc 
changes in a way that is not dependent on the presence of 
wild-type p53. 

In order to analyze the relationship between our metabohc 
analysis and transcriptional profiling, microarray datasets of 
MDA-MB-231 cehs and normal breast tissues were recovered 
from publicly accessible gene expression profile dataset (GEO 
and CellMiner). Data were analyzed as previously described 
(Balestrieri et aU 2009), and results indicated that several 
glycolytic genes (i.e., Hkl, Pfkp, Enol, and Ldha) were 
expressed at a higher level in cancer cehs as compared with 
normal breast tissue (Supplementary Figure S8A], in agree- 
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ment with the more active glycolysis shown in Figure 3B. 
Differences in gene expression were also observed for several 
TCA cycle genes (Supplementary Figure S8B} . 

These findings, collectively, suggest that oncogenic K-Ras 
transformation induces cell-specific transcriptional programs 
that have a role in the metabohc alterations we observed (i.e., 
increased glycolytic activity) and highlight the utility of 



systems-level approaches. Figure 4 depicts a schematic view 
of metabohsm, indicating the changes we observed upon 
K-Ras transformation using ^^C MFA and transcriptional 
profiling. In particular, it can be observed that both K-Ras 
transformed mouse and human cells, as shown by previous 
metabohc analyses, divert most part of consumed Glc to Lac 
(Figure 4A, red arrows) . Such a preferential metabolic flux is 
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supported by the observed high level of expression of 
significant glycolytic genes in both cell lines (i.e., Adpgk, 
Gadphy Pgky Eno 1 and 2, Pkm2 and Ldha; Figure 4 A, right 
box). In addition, the reduced flux of pyruvate entering into 
mitochondria (Figure 4 A, green arrow) and the deregulation of 
several genes encoding TCA cycle enzymes (Figure 4 A) in part 
provide an explanation of the decreased TCA cycle flux 
observed in both transformed cells lines (Figure 4A, green 
arrows). At the same time, our metabohc data indicate 
that such a reduced TCA cycle flux in transformed cells has 
an important anabohc role and is essentially fueled by Gin, as 
suggested by the identification of labeled Asp, a product of 
Gin-derived TCA cycle intermediates (Figure 4A). Further- 
more, the important role of Gin in supporting anabohc 
processes of cancer cells (Figure 4B, red arrows) is suggested 
by the identification of labeled metabolites containing 
Gin-derived nitrogen and associated increases in the expres- 
sion of genes involved in glutathione metabohsm [Gclm, Gss, 
and Ggtoi), hexosamine metabohsm (Gfptl), and nucleotide 
biosynthesis [Cad, Adss, and Adsl) (Figure 4B). Taken together, 
these findings indicate a major alteration in cancer cell 
metabohsm with a decoupling between Glc and Gin fates. 



Metabolic alterations are dependent upon 
oncogenic K-Ras activation and glutamine 
utilization 

To directly assess the role of oncogenic K-Ras in the previously 
described metabohc alterations of transformed cehs, 
we performed several metabohc analyses by taking advantage 
of a ceh line, NIH3T3-GEF-DN (reverted cehs, R), stably 
expressing a dominant-negative Guanine Exchange Factor 
CDC25 (GEF-DN) that specifically attenuates oncogenic K-Ras 
activation (Vanoni et al 1999; Bossu et aU 2000). Notably, R 
cehs exhibit changes consistent with a regression of cellular 
transformation, including Ras-GTP levels, morphology, 
anchorage-independent growth, reduction of Ras-dependent 
tumor formation in nude mice, glucose and glutamine 
dependence and mitochondrial dysfunction to a phenotype 
similar to that observed in N cehs (Chiaradonna et al, 2006b; 
Gaglio et al 2009). 

To compare the metabohc phenotype of R cehs with N and T 
cells, the intracellular concentration of Cit and Mai (Figure 5 A 
and B) or Glu and Asp (Figure 5C and D) was measured by 
enzymatic assays in cells cultured for 54 h in optimal growth 
conditions (25 mM Glc + 4mM Gin). Results indicated that 
N and R cehs contained similar abundances of these 
metabolites, while T-ceh metabohc abundances were similar 
to those previously observed using GC/MS. To further detail 
the metabohc characteristics of R cehs as compared with N and 
T cell lines and evaluate whether Gin utilization is a general 
requirement for the prohferation of K-Ras transformed cehs, 
we performed ceh prohferation analyses in medium containing 
low Gln+ each of two ceh permeable glutamine products, 
dimethyl aspartate (DMD) and dimethyl a-ketoglutarate 
(DMK). As shown in Supplementary Figure S9, Gin shortage 
caused a prohferation reduction only in T cehs, which was 
partially rescued by addition of DMD or DMK (Supplementary 
Figure S9). 
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To further test whether Gin metabolism induced in T cells was 
dependent on the K-ms oncogene, we treated N, T, and 
R cell lines with amino oxyacet ate (AOA), an inhibitor of 
aminotransferase activity (Figure 5E) and epigallocatechin 
gallate (EGCG), an inhibitor of glutamate dehydrogenase 
(GDH) activity (Figure 5F). As shown in Figure 5E, only T-cell 
prohferation was inhibited by AOA treatment (Figure 5E, 
compare middle panel versus left panel :N cells and right panel: 
R cells) and rescued by DMD addition and not by DMK addition 
(Figure 5E) . These data suggest that aminotransferase activity is 
necessary to sustain cancer cell growth and Gin-derived Asp is 
an important anabolic precursor for cancer cell growth. 

Analogous prohferation experiments were performed in N, 
T, and R ceh lines upon treatment with EGCG. As shown in 
Figure 5F, the treatment induced cell death in all cell lines. 
Nevertheless, the reduction of ceh number was markedly 
greater in T cells compared with that observed in N and R cehs. 
Importantly, addition of DMK completely rescued viability in 
ah the three ceh lines (Figure 5F) . Given the significant level of 
ceh death observed in ah cehs at 24 h, we investigated 
the effects of EGCG on ceh survival at earher time points 
(between 0 and 8 h) to better understand the mechanism of 
cytotoxicity. EGCG reduced viability in all three cell lines 
in a time-dependent manner, with T cehs being the most 
sensitive to EGCG treatment (Figure 6A, hght gray bars). 
Interestingly, the antioxidant N-acetyl-L-cysteine (NAG) 
protected only T cehs from EGCG-mediated cytotoxicity 
(Figure 6A), suggesting a specific role of ROS accumulation 
for their ceh death. To ascertain whether EGCG induced a 
different type of cell death in T cells as compared with N and R 
cehs, cehs were treated for 8h with EGCG and stained with 
Annexin V (AV) and propidium iodide (PI), indicators of late 
apoptosis and of plasma membrane permeability, respectively. 
Untreated cells exhibited an AV-negative/PI-negative pattern 
showing viable cehs (Figure 6B, upper panels), whereas N and 
R cehs treated with EGCG for 8 h exhibited an AV-negative/PI- 
positive pattern, showing a dead population with loss of 
plasma membrane integrity, typical of necrosis (Figure 6B, 
bottom panels). T cehs treated with EGCG contained a 
population of AV-positive/PI-positive pattern, showing a 
typically apoptotic population with an increased AV signal. 
These results were confirmed also by ceh morphology analysis 
after EGCG treatment. While EGCG induced flattened and 
attached ceh remnants in N and R cehs (almost ah Pl-positive 
cehs), EGCG-treated Tcehs were detached and rounded cehs 
(data not shown) . 

Given the fact that ATP depletion is more likely to induce 
necrosis, while ROS accumulation is more typically associated 
with apoptosis, we measured intracellular ATP and ROS levels 
upon 3 h of EGCG treatment in each ceh line. As shown in 
Figure 6C, EGCG caused a significant decrease of intracellular 
ATP levels in N and R cehs, which could not be rescued by NAG 
treatment. On the contrary, ROS levels were high and almost 
constant in T cehs regardless of EGCG treatment, whereas N 
and R cehs experienced a decrease in ROS levels upon EGCG 
treatment (Figure 6D) . Absolute levels of intracellular ROS in 
untreated Tcehs at 3 h were about two-fold higher as compared 
with N and R cells (e.g., N=193.5; 7=317.2; R=U1.7), and 
addition of NAG reduced ROS levels in ah cehs (Figure 6D). 
Altogether, these findings suggest that GDH activity, through 
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Figure 5 Relative metabolites concentrations in N, T, and R cell lines and effect of aminooxyacetate (AOA) and epigallocatechin gallate (EGCG) treatments on their 
proliferation. Evaluation of Cit (A), Mai (B), Glu (C), and Asp (D) intracellular concentrations was carried out after 54 h of grown in normal medium by using enzymatic 
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complete medium replacement, were treated with 100|iM AOA, or 1 mM dimethyl aspartate (DMD), or 2mM dimethyl a-ketoglutarate (DMK) or 20|iM EGCG and 
counted at indicated time points. Error bars indicate s.e.m. (n=3). 



different mechanisms is necessary for the survival of all three 
cell lines, and its inhibition may cause in N and R cells 
an energy derangement followed by necrosis, due to ATP 
depletion, and in T cells an antioxidant depletion that 
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ultimately may bring to apoptosis. In conclusion, these results 
suggest a role of glutamine in endogenous antioxidant 
production, for instance glutathione synthesis, further 
confirming the anabohc role of glutamine in T cells. 
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Figure 6 (A) Effects of EGCG inhibitor at sliort term in N, T, and R cell lines. Cells, seeded at density of 3000 cells/cm^, after 1 8 h and following a complete medium 
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analyzed after 8 h. The images were acquired by x 60 objective using Metamorph 7 software. (C) Analysis of intracellular ATP level by enzymatic assays. The cells 
were treated as above described and collected at indicated time points. (D) Intracellular ROS levels were measured by using 5 mM DCFH2-DA staining. The cells were 
treated as above described and collected at indicated time points. Error bars indicate s.e.m. (n=3). 



Discussion 

Cancer is a complex disease that arises from numerous 
molecular changes on a various cellular pathways (Mazurek 
et aU 1998; Hanahan and Weinberg, 2000). Among these 
alterations is the reprogramming of metabolic pathways 
associated with bioenergetics and cellular biosyntheses, and 
metabohc interventions are now emerging as potential 
therapeutic targets (Michelakis et aU 2010; Vander Heiden 
et aU 2010; Wise and Thompson, 2010). In this paper, we used 
^^C MFA, NTFD analysis, and transcriptional profiling to 
provide a comprehensive analysis of cancer cell metabohsm 
and of its regulation by the K-ras oncogene. The results, 
summarized in Figure 7, highlight the reprogramming of 
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metabohc pathways that occurs upon introduction of mutant 
forms of K-Ras. In cancer cells, most consumed Glc is diverted 
to Lac (Figure 7, red arrows), while the oxidation of Pyr to 
AcCoA via PDH is consistently decreased (Figure 7, green 
arrows) . This decrease results in a marked reduction of TCA 
cycle flux (Figure 7, green circle arrow), presumably caused by 
expression of Pdk genes. It is intriguing to note that both 
mouse and human transformed cells herein analyzed, present 
a substantially lower mitochondrial Complex I activity 
(Ishikawa et aU 2008; Baracca et aU 2010). Consistent with 
such a decrease in glucose utilization through the TCA cycle, in 
cancer cells Gin becomes the more relevant source of both 
carbon and nitrogen for cellular biosynthesis (Figure 7, black, 
blue, and sky-blue arrows), through aminotransferase 
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Figure 7 Schematic representation of Glc and Gin pathways decoupling in K-Ras cancer cells. Glc provides the main energy source for cancer cells with oncogenic K- 
Ras. Instead, Gin provides the carbon and nitrogen for the synthesis of biomass building blocks (amino acids, nucleotides, and fatty acids). Red arrows represent 
enhanced flux, green arrows represent reduced flux, black arrows represent no change flux, blue arrow represents fonward flux of TCA cycle, and sky-blue represents 
reductive carboxylation flux. 



(Asp labeling) and glutamate dehydrogenase activities 
(Figure 7, black, blue, and sky-blue arrows). 

In conclusion, in K-Ras transformed cells the decoupling 
between Glc and Gin metabohsm leads to an efficient 
utilization of both carbon and nitrogen from Gin into biomass 
building blocks (amino acids and nucleotides) and glu- 
tathione. Both, sustaining growth and the ability to quench 
ROS production, may promote cancer cell prohferation and/or 
survival. The strong dependence of transformed cell growth on 
Gin availability is, therefore, not surprising (Gagho et aU 2009; 
Wang et aU 2010; Weinberg et aU 2010). Our analysis also 
suggests that transcriptional changes contribute to the 
observed metabohc reprogramming in mouse and human 
transformed cells. In the future, detailed analyses of the 
various enzymatic steps involved in the growth-sustaining 
metaboHsm as well as genes regulating these reactions 
may shed light on the molecular basis of transformation. In 
addition, these results may enable the identification of fragile 
points in the metabohc network and more effective, targeted 
cancer therapeutics. 

Materials and methods 
Cell culture 

Mouse embryonic fibroblast NIH3T3 cells (CRL-1658; American Type 
Culture Collection) and a K-Ras transformed NIH3T3 -derived cell line, 
226.4.1 (Pulciani et at, 1985), were routinely grown in Dulbecco's 
modified Eagle's medium containing 10 % newborn calf serum, 4 mM 
L-glutamine, lOOU/ml penicillin, and lOOmg/ml streptomycin (nor- 
mal growth medium) (all from Invitrogen, Carlsbad, CA, USA) at 37°C 
in a humidified atmosphere of 5 % CO2. The NIH-GEF-DN (reverted cell 
line), stably and constitutively expressing the dominant-negative 
mutant Cdc25MmW1056E (GEF-DN; Bossu et at, 2000), was main- 
tained in normal growth medium supplemented with 0.7mg/ml 
geneticin (G418; Sigma- Aldrich Inc., St Louis, MO, USA) . MDA-MB-231 
human breast carcinoma cell line, obtained from the European 
Institute of Oncology (Milano, Italy), was grown in Dulbecco's 
modified Eagle's medium containing 5% fetal bovine serum, 4mM 
L-glutamine, lOOU/ml penicillin, and lOOmg/ml streptomycin 
(normal growth medium) (all from Invitrogen, Carlsbad, CA, USA) 



at 37°C in a humidified atmosphere of 5% CO2. HCT116 p53 + /+ and 
HCT116 p53~^~ cells, obtained from the European Institute of 
Oncology (Milano, Italy), were cultured with complete DMEM plus 
10% FBS. Cells were passed using trypsin-EDTA (Invitrogen) and 
maintained in culture for 48 h before manipulation. 



Cell proliferation analysis and cell treatments 

For long-term studies, cells were plated at the density of 3000 cells/cm^ 
in complete growth medium. After 18 h, the cells were washed twice 
with phosphate-buffered saline (PBS) and incubated in a media with 
different glucose concentrations (25 and 1 mM glucose) or different 
glutamine concentrations (4 and 0.5 mM glutamine) (time 0). 
When indicated, cells were treated with 100 |iM AOA (0-(Carbox- 
ymethyl)hydroxylamine hemihydrochloride), or 20|iM (-)-EGCG, or 
2 mM dimethyl 2-Oxoglutarate -DMK- or 1 mM L-Aspartic acid 
P-methyl ester hydrochloride -DMD- (Sigma- Aldrich Inc.) . For short- 
term studies, cells were plated and grown as above and then treated 
with 20 |iM EGCG or 5 mM NAC (Sigma-Aldrich Inc.). To measure cell 
proliferation, harvested cells were counted by Burker chamber. 



Microarray data collection 

Microarray transcriptional profiles were extracted from two different 
datasets. The mouse time course data were obtained from a 
transcription profiles collection of NIH3T3 and NIH3T3 K-Ras mouse 
fibroblasts, generated in our laboratory, grown in 25 mM gluco- 
se + 4 mM glutamine, along a time course of 72 h. In particular, the 
cells used for the time course transcription analysis, were collected 
at 18h from the initial seeding, time corresponding to the change of 
medium (indicated as To) and at 24, 48 and 72 h (NCBI GEO database 
accession from GSM741354 to GSM741361 and from GSM741368 to 
GSM741375 of NIH3T3 and NIH3T3 K-Ras, respectively) . For each time 
points, labeled cRNA was generated by using the Affymetrix One-Cycle 
Target Labeling and Control Reagent kit (Affymetrix Inc., Santa Clara, 
CA, USA), following the manufacturer's protocol and was hybridized 
using Affymetrix Genechips (Mouse Genome 430 2.0 Array) in order to 
determine the global gene expression patterns. Arrays were washed 
and scanned on the Affymetrix Complete GeneChip® Instrument 
System and processed into GEL files. 

Human raw expression data for U133A arrays were obtained from 
publicly accessible gene expression profile database GEO (http:// 
www.ncbi.nlm.nih.gov/geo/) or CellMiner (http://discover.nci.nih. 
gov/cellminer/home.do). In particular, we used the transcription 
profiles data of the human breast cancer cell line MDA-MB-231 
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(CellMiner accession 13409hgl33a21) and of the human normal breast 
tissue (NCBI GEO accession GSM44683). Both mouse and human 
transcriptional data were separately submitted on the same procedure 
of normalization and filtering of noise, in order to obtain a good degree 
of comparison across the different experiments. Briefly, the data files 
(CEL files) were imported into GeneSpring GX 11.0.2 software (Agilent 
Technologies Inc.), by using RMA method, were normalized and 
summarized as probe-level intensity measurements (Irizarry et at, 
2003). In particular, the algorithm implemented in GeneSpring, starting 
from the perfect-match probe-level data of a set of arrays, performs the 
background correction, the normalization and finally summarize the 
results as a set of expression measures for each probe set. Subsequently, 
also 'per-gene normalization' was performed as described in Gene- 
Spring's manual. Altogether, these procedures permitted to obtain the 
relative expression of each probe set in a log2 scale with the value 0 as 
center. This scale of measurement was used for the heat maps showed 
in the paper. From these two datasets (mouse and human) have been 
identified and gathered the expression levels for genes encoding 
proteins involved in glycolysis, TCA cycle and glutamine metabolism, 
at 48 and 72 h for mouse datasets and, at time indicated by the authors 
for both human tissue and cell line samples. Total selected mouse and 
human genes and their intensity expression levels have been reported 
in Supplementary Dataset. 



Metabolite extraction 

For labeling experiments and GC/MS-based metabolomics, asynchro- 
nous NIH3T3, NIH3T3 K-Ras, and MDA-MB-231 were plated at 3000 
cells/cm^ in 6-well plates with normal growth medium. In all, 18 h 
after seeding, cells were washed with PBS and incubated in different 
DMEM (Sigma) plus 10% dialyzed serum supplemented with 
unlabeled glucose or glutamine and either 4mM [U-^^Cs] glutamine 
or a 25 mM 1:1 mixture of [U-^^Cs] glucose and [l-^^C]glucose. After 
54 h, spent medium was collected and analyzed for glucose, lactate, 
and glutamine consumption on a YSI7100 analyzer. Cells were quickly 
rinsed with PBS to minimize spent medium carryover and quenched 
with 0.4 ml ice-cold methanol. An equal volume of water containing 
1 |ig norvaline internal standard was added, and cells were collected by 
scraping with a pipette tip. Two volumes of chloroform were added, 
and the cells were vortexed at 4°C for 30min. Samples were 
centrifuged at 3000 g for 10 min, and the aqueous phase was collected 
in a new tube and evaporated under airflow at room temperature. 



Derivatization and GC/IVIS measurements 

Dried polar metabolites were dissolved in 60 |il of 2 % methoxyamine 
hydrochloride in pyridine (Pierce), sonicated for 30 min, and held at 
37°C for 2h. After dissolution and reaction, 90|il MBTSTFA+1% 
TBDMCS (Pierce) was added and samples were incubated at 55°C for 
60 min. GC/MS analysis was performed using an Agilent 6890 GC 
equipped with a 30-m DB-35MS capillary column connected to an 
Agilent 5975BMS operating under electron impact (EI) ionization at 
70 eV, 1 |il of sample was injected in splitless mode at 270°C, using 
helium as the carrier gas at a flow rate of 1 ml/min. The GC oven 
temperature was held at 100°C for 3 min and increased to 300°C at 
3.5°C/min for a total run time of ~ 60 min. The MS source and 
quadrupole were held at 230 and 150°C, respectively, and the detector 
was operated in selected ion monitoring mode. MIDs were obtained for 
each measured metabolite and incorporated with extracellular flux 
measurements for flux determination. For determination of relative 
metabolite abundances, the integrated signal of all potentially labeled 
ions for each metabolite fragment was normalized by the signal from 
norvaline and the per well cell number (obtained by counting 
surrogate plates) . 



Flux estimation 

Intracellular fluxes were estimated for a model reaction network by 
minimizing the lack of fit between actual and simulated flux and 
GC/MS measurements. The network contained simplified versions of 
glycolysis, the PPP, anaplerotic reactions, the TCA cycle, and amino- 



acid biosynthesis. Specific assumptions regarding the network and 
metabolism that were made for the sake of the flux analysis are listed 
in Supplementary information. We additionally calculated 95% 
confidence intervals for each flux using parameter continuation 
(Antoniewicz et at, 2006). All flux simulation, estimation, and 
continuation in this study were conducted using Metran, a flux 
analysis tool built upon an EMU framework (Antoniewicz et at, 2007; 
Yoo et al, 2008; Young et al, 2008) . When required, MIDs were corrected 
for natural isotope abundance using in house algorithms adapted from 
Fernandez et al (1996). Simulated and measured MIDs for MFA of N and 
T cells are presented in Supplementary Figures S2 and S3 ([U-^^Cs] 
glutamine) or SIO and Sll ([U-^^C6]glucose+ [l-^^C]glucose). Data 
for MDA-MB-231 cell MFA ([U-^^Cs] glutamine) are presented in 
Supplementary Figure S6. 



Non-targeted tracer fate detection 

The NTFD analysis was performed to determine the metabolic fate of 
the a-amine nitrogen atom of glutamine (Hiller et al, 2010). Cells were 
incubated for 54 h as described above, but instead of unlabeled 
glutamine, a 1:1 mixture of unlabeled and [ a- ^^N] glutamine (Cam- 
bridge Isotope Laboratories) was employed for each growth condition. 
Since it is a requirement of NTFD, additional cells were cultivated with 
completely unlabeled glutamine as reference. Both cultivations were 
performed in triplicate and intracellular metabolites were extracted 
and measured as described above. The raw GC/MS data were used as 
input for the NTFD software and all detected labeled compounds for all 
growth conditions were matched using algorithms of the Metabolite 
Detector software package (Hiller et al, 2009). Finally, the NIST 
reference library was employed for compound identification (Babush- 
ok et al, 2007). 



Metabolite assay 

For enzymatic assay quantitation of aspartate, citrate, malate, and 
glutamate, asynchronous NIH3T3, NIH3T3 K-Ras, and NIH-GEF-DN 
were plated at 3000 cells/cm^ with normal growth medium. In all, 18 h 
after seeding, cells were washed with PBS and incubated in medium 
with 25 mM glucose and 4mM glutamine. After 54 h, the cells were 
washed twice with PBS cold and collected using trypsin-EDTA 
(Invitrogen) . The cells were counted and 500 000 cell aliquots were 
rapidly homogenized with appropriate assay buffer on ice. The 
samples were centrifuged at 15 000 g for 10 min to remove cell debris 
and deproteinized using a perchloric acid/KOH protocol (Bio Vision) . 
The samples were added into triplicate wells of a 96-well plate and the 
volume was brought to 50 |il with assay buffer. The samples were 
incubated for 30 min at room temperature or at 37°C, protected from 
light. Then, OD was measured in a microplate reader. 



Intracellular ATP quantification 

Intracellular ATP levels were measured using CellTiter Glo® luciferin- 
luciferase assay (Promega) according to the manufacturer's protocol. 
In particular, 100 [i\ of cell suspension containing 10^ cells was added to 
an equal volume of CellTiter Glo Reagent in a single well of 96-well 
plate. The emitted luminescence was collected after incubation at the 
wavelength of 560 nm using Gary Eclipse spectrofluorimeter (Varian) 
and the luminescence values were converted in ATP quantities after 
the setting of a calibration curve. 



ROS levels measurement 

Detection of ROS levels was carried out through flow cytometric 
analyses using a FACScan flow cytometer (Becton-Dickinson) with 
CellQuest software (BD Biosciences). ROS levels were measured 
through staining plated cells with 5 |iM dichloro-dihydro-fluoresceine- 
diacetate (DCFH2-DA; Molecular Probes, Invitrogen) for 30 min at 
37°C. Then, the cells were trypsinized, collected in ice and acquired by 
FACScan. Analysis of flow cytometric data was carried out using 
WinMDI software. 
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AV/PI cell death assay 

Cell death was estimated using the Apoptosis Detection Kit (Immunolo- 
gical Sciences). Cells were plated on coverslips previously treated with 
2 % gelatin. Cells on coverslips were stained by adding 50 |il of binding 
buffer, 5 |il of Annexin V-FITC, and 5 [xl of PI. After 15 min of incubation at 
37° C in the dark, cells were treated with 5 |ig/ml Hoechst (Invitrogen) for 
5 min at RT. The cover glasses, mounted in DABCO, were analyzed under 
a Nikon ECLIPSE 90i fluorescence microscope equipped with a b/w CCD 
camera (Hamamatsu-CoolSNAP, Hamamatsu Corporation Japan), using 
Plan Apo objective ( x 60 oil; numerical aperture 0.75 and 1.4, 
respectively). The images were acquired using the imaging software 
Metamorph 7, and then processed in Adobe Photoshop CS3 with 
adjustments of brightness and contrast. 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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